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Abstract. In multi-level systems, the commonly used adiabatic elimination is a 
i— i method for approximating the dynamics of the system by eliminating irrelevant, non- 

resonantly coupled levels. This procedure is, however, somewhat ambiguous and it is 
I 1 not clear how to improve on it systematically. We use an integro-differential equation 

for the probability amplitudes of the levels of interest, which is equivalent to the original 
Schrodingcr equation for all probability amplitudes. In conjunction with a Markov 
approximation, the integro-differential equation is then used to generate a hierarchy of 
approximations, in which the zeroth order is the adiabatic-elimination approximation. 
It works well with a proper choice of interaction picture; the procedure suggests criteria 
for optimizing this choice. The first-order approximation in the hierarchy is found to 
be sufficient for practical purposes, and is not so sensitive to the choice of interaction 
^sO picture. We illustrate the method for a single three-level atom and a pair of three-level 

atoms with Rydberg blockade. 
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1. Introduction 



The coherent manipulation of quantum states is central to a large variety of contempo- 
rary research in physics, among them experiments that aim at processing quantum 
information. When implementing a qubit in the internal degrees of freedom of a 
trapped atom or ion, the relevant quantum states are different electronic energy levels 
manipulated by a combination of static and time-dependent electromagnetic fields. The 
ideal case of a two-level system resonantly driven by a laser is rarely available as the levels 
of interest may not be coupled directly by a dipole transition, or the transition occurs 
at an inconvenient wavelength. Instead, two levels are often coupled indirectly through 
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intermediate levels, leading to a multi-level system, which can be very complicated to 
analyze exactly. Fortunately, the intermediate levels are usually barely populated and 
are of little consequence, other than providing a coupling route, to the dynamics of the 
levels of interest. It is then common to employ the procedure of adiabatic elimination 
to remove these irrelevant levels, thereby achieving a considerable simplification of the 
formalism in terms of a lower-dimensional effective description. 

In this paper we first review the standard method of adiabatic elimination and the 
problems that may arise when applying it (Section [2|. Then, in Section [3j we present a 
systematic approach for separating the more relevant and the less relevant subspaces by 
means of an integro-differential equation of Lippmann-Schwinger type that can be used 
to generate a hierarchy of approximations, starting with the usual adiabatic elimination. 
A technical point of considerable importance is the choice of an appropriate interaction 
picture; we propose optimality criteria for a judicious choice. In Section [4] finally, we 
use a driven Raman transition of A-type in a three-level atom, and also two atoms with 
a three-level cascade experiencing a Rydberg blockade, for case studies that illustrate 
the usefulness of our approach. 

2. Brief review of adiabatic elimination 

2.1. Example: Raman transition in a three-level system 

In a three-level system, stimulated Raman transition is a two-photon process used to 
couple a ground state |0) to a target state |1) via a far-detuned intermediate state |e), 
with respective energies huo,hu!i, and hu e . The use of such an indirect coupling is 
necessary whenever a direct transition between the ground state and the target state 
is not possible, because the transition is dipole forbidden, or has a transition frequency 
outside the range accessible to common lasers. Depending on the energies of the three 
levels, the system is either in the A configuration (u ,Ui < u c ), the V configuration 
(u>o,u>i > co e ), or the cascade configuration (wo < u c < ui). 

To be specific, we will focus on the A configuration of figure [TJ but the other 
configurations can be dealt with in very much the same way. The transitions are driven 
by two distinguishable lasers (with different frequencies or polarizations, for example) 
of frequencies o;lo an d Wli- The ith laser has a detuning Aj = u c — oJi — uu to the 
intermediate state. It is convenient to introduce the average detuning A = ^(A + Ai), 
as well as the overall detuning 5 = A — Ai of the two-photon process. 

After applying the dipole and rotating-wave approximations, the Hamiltonian of 
the three-level system is represented in the basis {|0), |1), |e)} by the matrix 

/ |n e^ L »*\ 

H = h wi |Oie iWL1 * , (1) 

y±tt*e^ ot ^*e~ iWL1 * u c i 

where the Rabi frequencies depend on the amplitude and the polarization of the 
electric field, as well as the dipole matrix elements between the states. In the interaction 
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Figure 1. The average detuning A is large compared with the Rabi frequencies of the 
lasers with frequencies wlo and wli that drive the transitions |0) — » |1) and |0) — » |1), 
respectively. The overall detuning S of the two-photon transition is small. 



picture denned by the splitting 



H = H + H 1 with H = h 
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the resulting Hamiltonian is time-independent 
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where we indicate the splitting of the 3x3 matrix into the top-left 2x2 matrix for the 
relevant levels |0) and |1) and the rest that involves the intermediate level |e). 

The time evolution of such a three-state problem can, of course, be studied 
using an exact eigen-decomposition, but the actual expressions for the eigenvalues 
and eigenvectors of Hi are often involved and not transparent. Therefore, adiabatic 
elimination is frequently used to remove the intermediate state from the description 
and thus reduce the three-level system to a two-level system of the relevant ground and 
target states, with a 2 x 2 matrix for its effective Hamiltonian. Then all the familiar 
tools for two- level systems are applicable, and one can easily work out the frequency and 
amplitude of the oscillatory probability amplitudes as well as other details of interest. 

Adiabatic elimination relies on the fact that the intermediate level, which is not 
populated initially, is only weakly and non-resonantly coupled to the ground state and 
the target state, i.e., |Oj| <C |A|. It is then customary to set, in the equations of motion 
in the interaction picture, 



i^Mt) = ln* c (t) + l&Mt) + Ac c (t) « 0, 



(4) 



<9T CVV 2 "" u ~ uw ' 2 
where c , cr, and c e denote the interaction-picture probability amplitudes for the respec- 
tive three levels depicted in figure [TJ Under this approximation, c e (t) can be expressed 
as a linear combination of Cq(£) and Ci(t), which then permits the elimination of c e (t) 
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from the equations of motion for the probability amplitudes Co(i) and ci(t). This leads 
to an effective two-dimensional Hamiltonian 



now written in the basis {|0), |1)} involving the relevant states only. 
2.2. Problems with adiabatic elimination 

Despite its popular use in simplifying multi-level problems, the procedure of adiabatic 
elimination has certain ambiguities and potential problems. We present these in the 
form of four questions: 

(a) The basic assumption that the time derivative of c e (t) is small and can be 
treated as approximately vanishing is questionable. Even if the population in the 
intermediate state is small, its time derivative need not be small; in fact, it is rather 
sizeable because the large detuning makes c e (t) oscillate rapidly. Is there a better 
justification for this assumption ? 

(b) There are problems with the normalization of the wave function. Originally, we 
have |co(t)| 2 + |ci(t)| 2 + |c e (t)| 2 = 1. As if e ff is hermitian, the two-component wave 
function that goes with it has |co(t)| 2 + |ci(t)| 2 = constant. Combined with the 
initial condition (co,ci,c e ) = (1,0,0) this implies c e (t) = for all t, whereas the 
basic approximation in ^ says otherwise. Thus it fails to make a consistent 
statement about the population in the intermediate state. Is it possible to 
estimate the population in the eliminated state and comply with the normalization 
conditions? 

(c) It is always possible to go to another interaction picture by subtracting 
a multiple of the identity from Hq, and split the Hamiltonian differently, 
H = (H — E) + (Hi + E) = H + Hi. This yields a different interaction-picture 
Hamiltonian, Hi = Hi + E. The modification appears innocuous, but it results in 
a different effective Hamiltonian upon the adiabatic elimination of the intermediate 
level. There have been attempts to identify the "correct" choice of interaction 
picture, which gave valuable insights pQ, but a definite answer is still lacking. How 
can one choose an optimal interaction picture? 

(d) For larger Rabi frequencies or smaller detunings, one observes that the adiabatic 
elimination does not give a trustworthy approximation. Is there a systematic way 
of improving the accuracy of the adiabatic- elimination approximation? 

We offer answers to these four questions in the next sections. 




(5) 
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3. General method 

3.1. Integro- differential equation for the relevant states 

In the more general situation of a multi-level system with m + n energy levels, we 
distinguish the set of m relevant states, whose probability amplitudes are collected in 
the m-component column ip(t), and the set of n auxiliary states, for which we have the 
n-component column e(t), with weak or highly non-resonant coupling between states 
from different sets. In the spirit of adiabatic elimination, we wish to have an effective 
description for the relevant states described by ip(t). 

It is most convenient to work in an interaction picture where the Hamiltonian is 
time- independent. For a basis where the first m elements correspond to the relevant 
states and the following n elements to the auxiliary states, the interaction-picture 
Hamiltonian is composed of the submatrices u, Q, and A, 



« = » w K l *(«)= IX (6) 





is the (m + n)-component column of probability amplitudes. The Schrodinger equation 
then takes the form 

iJU(» =e*J>(t) + iW), 

Ot Z 

= inty(f) + Ae(f). 

de 

If we were using adiabatic elimination, we would now set — = 0. Instead, we solve the 
differential equation for e(t) with the initial condition of no population in the auxiliary 

states, e(t = 0) = 0, 

t 

e(t) = ~\J dt'e~ iA ^n^(t'). (8) 
o 

We then use this in the differential equation for ip(t) to arrive at an integro-differential 
equation of Lippmann-Schwinger type, 

t 

iJ^W = ^(t) dt'e'^-^nmt'). (9) 

o 

Together with (|8p, this is fully equivalent to the pair of equations in (|7|. While no 
approximation entered in the transition from ^ to ([9]), the integro-differential equation 
now serves as the starting point for the generation of a hierarchy of approximations. 
This approach is similar in spirit to, but differs in detail from, the Green's-function 
approach discussed in pp. 
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3.2. Markov approximation 

The usual (zeroth-order) Markov approximation assumes that memory effects are 
negligible, so that ip(t') ~ ip(t) in the integral in ([8]). This assumption is valid, provided 
that tp{t) oscillates much more slowly than e _lA *. This happens when |0|, \lo\ <C |A|, 
which is the case when the auxiliary states are coupled weakly and non-resonantly to the 
relevant states. As discussed in standard texts, such as references [2] and [3], a coarse 
graining is then appropriate in the evaluation of the remaining integral over t' because 
time scales of order A -1 are not resolved and e _lA * averages out over the coarse-grained 
time intervals. Accordingly, we have 
t 

df / e -iA( t -f) = — w _L (io) 

iA iA 

and thus obtain the effective Hamiltonian for the zeroth-order Markov approximation, 

H<g = hL-n±rt\ (ii) 

This happens to be exactly the effective Hamiltonian that adiabatic elimination would 

yield. In this sense, then, we have given an answer to question daj) in Section 2.2 Setting 

d — 

— e e (i) ~ in the adiabatic elimination procedure amounts to a shorthand for the above 

coarse graining in time to discard rapidly oscillating features and retain only the slowly 
varying dynamics. 

To gain some information about the evolution of e(t), and so arrive at an answer to 
question |b]), we consider the first-order Markov approximation that takes a bit of the 
history of the relevant states into account in the integral in (|8]). We accomplish this by 
means of the approximation 

m^m-it-f) 9 ^, (12) 

from a Taylor-series expansion about t' = t. With (10) and 

t t 

dt' (t - i')e- iA ^ = J dt'e- iA ^ « (13) 
o o 
this amounts to 

dt V 4A / 4A 2 dt K J 

In defining the effective Hamiltonian, we now have a choice between 

(l + a^y^-a^m^HSlm (is) 

and 

i»f> = (i + tij^m) " 2 («, - (i + n ^ n .) "'Vw 

= jrfffiVM (16) 
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with 



1/2 



(17) 



The effective Hamiltonian in (16) describes the evolution of the scaled wave function 



ijj'(t), and is hermitian in the usual sense, so that 



(18) 



is constant in time. By contrast, the effective Hamiltonian in (15) is hermitian only for 
the modified inner product 

1 



(fc^) = ^i + fl-fitU 



4A 2 



-2: 



(19) 



giving the right-hand side of (18) for the normalization of ip(t) in accordance with 
(■?/>, = 1- It does not matter which formulation we prefer for the first-order Markov 
approximation, either (15) with (19) or (16) with (17) and (18) is fine. The two 



formulations are related to each other by the similarity transformation afforded by 
(1 + |fiA _2 f2t) 1 / 2 ; anc i reliable results can only be expected if this operator does not 
differ much from the identity. This gives a more precise meaning to the requirement 
that the coupling part Q, of the interaction-picture Hamiltonian in ^ should be small 
on the scale set by the energy part A of the auxiliary levels. 

A comparison of (18) with the normalization of the full column \I/ of (J6|), 

*t^ = ^ + e t e= i j (20) 
reveals that ee is here approximated by 

^n-^nty. (2i) 



e'e 



4A 2 



This is just what the zeroth-order Markov approximation says, 



e(f) 



2A 



(22) 



This observation is reassuring and provides the answer to question ^ in Section 2.2 



The first-order Markov approximation offers a correction to adiabatic elimination. 



By including another term in the Taylor-series approximation (12), one gets a second- 



order Markov approximation. A second time derivative of ip(t) will appear, but it can be 
approximated by replacing one of the time derivatives by an application of the effective 
Hamiltonian, thereby obtaining a Schrodinger-type first-order differential equation for 
ip(t). Successive terms in the Taylor series give a systematic hierarchy of approximations, 
with adiabatic elimination as the zeroth-order procedure, and this offers an answer to 
question (Jd|). As we shall see shortly, however, the first-order Markov approximation is 
already of sufficient quality, so that the complications of higher-order approximations 
are hardly worth the trouble. 
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3. 3. Choice of interaction picture 



There still remains question Q about the choice of a good, or perhaps optimal, 
interaction picture that leads to an effective Hamiltonian with the correct physics. Since 
we prefer interaction-picture Hamiltonians that are constant in time, Hamiltonians for 
different choices of interaction picture can only differ by a multiple of the identity. 



The effective Hamiltonians and H^' of (|11|) and (|15|), respectively, in an alternate 
interaction picture, would therefore take the 



orm 



H. 



(0) 



off 



hi u + u — fi 



4(A + u 



-fit 



H. 



(i) 



cff 



1 + Q- 



:fi f 



-1 



(23) 



H. 



(o) 

eff ' 



4(A + w) 

where to = E/h for the energy parameter E in question (|c]). 

We recall that the crucial assumption in the Markov approximation is that ip(t) 
oscillates slowly in comparison with e _lA *, which is invoked when replacing ^(t') by 
ip(t) on the right-hand side of (12) in the integro-differential equation Q. Consistency, 
therefore, requires that the magnitudes of the eigenvalues of the effective Hamiltonian 
should be as small as possible because they determine the frequencies contained in ip{t). 
The following three conditions for choosing the value of cu suggest themselves: 

(a) tr{u; + u} = 0, (b) ||i? e ff||op is minimal, (c) ||i? e ir||tr = tr | -H" c ff | is minimal. (24) 
Here, \\H || op refers to the operator norm of H, while || tr is the trace norm. Conditions 



(24 3) and (24 2) insist on choices of u such that if e ff has eigenvalues that are small in two 



different senses: The operator-norm condition requires minimizing the largest absolute 
value of the eigenvalues of H c g] the trace- norm condition requires minimizing the sum of 
the absolute values of the eigenvalues of H c g. Condition (24i) prescribes distributing the 
eigenvalues of the u> part of the interaction-picture Hamiltonian (JsT) about 0. This does 



not necessarily guarantee small eigenvalues, but for the Hamiltonians in (23), a different 



choice of u gives, as the dominant change, only an overall shift in the eigenvalues of 



the Hamiltonian. Condition (24i) is hence a simple alternative to a condition based 
on norms of H e ff, and occurs in the context of adiabatic elimination — the u of the 
interaction-picture Hamiltonian Hj in ^ that gave the effective Hamiltonian ^ from 



adiabatic elimination satisfies condition (24i). The example in Section 4.1 provides 



some evidence in favor of conditions (24d) and (24 2), with minimal difference between 



them, but the quality of the approximation is not substantially worse if condition (24i) 
is enforced. 



4. Approximation of Three-Level Systems 

4-1- One atom 

In order to check the validity of, and illustrate the differences between, the adiabatic- 
elimination approximation and the first-order Markov approximation, we apply them to 
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a three-level system of the kind depicted in figure [T} For a quantitative statement 
about the quality of the approximations, we compare the resulting Rabi frequency 
f^R for population transfer between |0) and |1). Here, the exact M7 R is the smallest 
spacing between the eigenvalues of Hi, and is approximated by the difference of the two 
eigenvalues of the H e s under consideration. We first consider the case of zero overall 
detuning. Expanding the solutions for 5 = to order x 4 , where x = |(|^o| 2 + |^i| 2 )/A 2 , 
gives 

' x 2 + 2x 3 - 5x 4 + 0(x 5 ) 



exact solution: 



n 



R 



Oth-order, (24i): — fi R 



Oth-order, (24d) 



Oth-order, (24 



A 
1 

A 
1 

A 
1 



x 



x. 



-(yr+Ax 



R 



X — 



1 2 , 1 . 
-X + -X 

2 2 



lst-order, (24i): — fi R 



lst-order, (124 



lst-order, ((24 



A 
1 

A 
1 

A' 



X 



X 



X 



x 2 + 2x A - 5x 



-x 4 + 0{x z 
+ 0(x 5 ) -- 



x 2 + x 3 



X 



+ 0(x 5 



1 

2 

x 



Vl + 2x - 1, 
(Vl + 4x-l). 



(25) 



1 + x 



7 o 15 
-x - 



-0 



R 



x — X 



4 4 
2 + x 3 - x 4 + 0(x 5 



0(x 5 



X 

1 + x' 



Regardless of the choice of conditions (245,)-(24^), the zeroth-order Markov 
approximation gives the correct Rabi frequency up to order x, while the first-order 
solutions give accuracy up to order x 2 . With a judicious choice of interaction picture, 
one can improve the approximations — for example, under the zeroth-order Markov 
approximation, employing condition (24:) in fact gives the exact Rabi frequency. 



However, we believe this is a coincidence rather than a general fact (see, for instance, 
the 5^0 case discussed below). A more reliable conclusion one can draw is that the 
first-order Markov approximation gives better accuracy, and that whether we choose 
condition (24i), (24d), or (24 z) matters little. In fact, one might even consider 



minimizing ||-ff e ff|| f° r other norms (for example, the Hilbert-Schmidt norm). The near- 
equivalence between the conditions is reassuring, since very often, condition (24i) is the 
easiest among the three possibilities to apply. It also validates the common choice of 
condition (24 a,) in standard adiabatic elimination. 

For comparison, we consider a case where 5^0. In many experiments targeting 
full population transfer between the ground state and the target state, 6 is adjusted to 
compensate for the light shifts of the atomic levels, due to the presence of the lasers, 
to give a resonant two-photon transition. Often, the choice is 5 = (|^i| 2 — |f2 | 2 )/4A, 
a value designed to make the diagonal entries of the effective Hamiltonian ^ from 
adiabatic elimination equal. This ensures a resonant process, to the accuracy justified 
by the adiabatic-elimination approximation. For this value of 5, the predicted Rabi 
frequencies for the exact solution and for the zeroth-order Markov approximation under 
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80 100 120 



Figure 2. Population probabilities for a single-atom Raman transition in a three-level 
system as a function of At, for Rabi frequencies Hq — 0.4A and Sli = 0.3A, and overall 
detuning S = (|^i| 2 — |Slo| 2 )/4A. The curves that start close to 1 show the population 
of the ground state |0); those that start at and increase to about 1 before decreasing 
again are for the target state |1); those that start near and never grow to large values 
are for the intermediate state |e). The solid black curves plot the exact solution; the 
thin dashed blue curves show the zeroth-order Markov approximation with condition 



(24i) (standard adiabatic elimination); the thick dashed blue curves also show the 



zeroth-order Markov approximation, but with condition (24:); the dotted red curves 



are for the first-order Markov approximation with condition ( 24 1) . 



the different conditions are 
1 



exact solution: 



Oth-order, d24 



Oth-order, (24d) 



Oth-order, p4 



A 
1 

A 
1 

A 
1 

A 1 



n 



R 



VT^{x - x 2 + [2 + 0(a 2 )]x 3 + 0(x 4 )}, 



n 



R 



a 2 x, 



a" 



x - -x 2 + 



+ 0{a 2 



x 3 + 0(x 4 



(26) 



R 



where we have defined the ratio a 



Vl -a 2 {x - [1 + 0(q 
|tto| 2 -|tti| 2 



2 ^x 2 + [2 + 0(or)\.r 



0. As in the 5 = example in (25) 



|ft | 2 + |fii| 2 

the zeroth-order approximation gives the correct answer to linear order in x, regardless 



of the choice of condition. That condition (24:) gives the exact Rabi frequency is no 



longer true here, although for small ft, it still provides the most accurate prediction. 
One can also verify that the first-order Markov approximation again gives the correct 
value for the Rabi frequency up to second order in x. A pertinent remark here is that, in 
this case where 5^0, conditions (|24b) and (24:) are rather more complicated to impose 



than condition (24i). This lends justification to improving the description of the system 
by going to higher-order Markov approximations using only the simplest condition (24 1), 
instead of correcting the zeroth-order approximation by a more complicated choice of 

UJ. 
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To further illustrate matters, we examine the predictions for the evolution of the 
populations in the three states under the different approximations for the case with Rabi 
frequencies flo = 0.4 A and Qi = 0.3 A, and overall detuning 5 = (|f2i| 2 — |fi | 2 )/4A. The 
results are plotted in figure [2j Standard adiabatic elimination (thin dashed blue lines in 



the plot), equivalent to the zeroth-order Markov approximation with condition (24i) 



gives a period for the Rabi oscillation that is about 6% short of the exact value. The 



zeroth-order Markov solution with condition (24 2), on the other hand, gives the exact 
Rabi frequency, but makes no statement about the population of the intermediate state. 
In contrast, the first-order Markov approximation, even with the simplest condition 



(24i), performs very well, and provides a reasonable estimate of the intermediate state 
population. Note that all the approximations show the effect of the coarse graining: 
They do not reproduce the high-frequency modulation of the exact solution. This coarse 
graining is also behind the fact that the population of the ground state, under the first- 
order Markov approximation, does not start at 1 at time t = 0. Instead, it begins with 
a value that can be thought of as the average over the initial coarse-grained time step, 
and is consistent with having an initial nonzero population in the intermediate state in 



compliance with (21) 



4-2. Two atoms 



For the single three- level atom in Section 4.1, one hardly needs an approximate 
treatment. Nevertheless, the reduction of the description to the relevant two levels 
is a helpful simplification and the coarse-grained probabilities thus obtained are often 
all one needs to know to determine the important parameter range for an experiment. 
When more levels are involved, a full treatment may no longer be possible or contain too 
many irrelevant details. As an example of such a more complex situation, we consider 
two identical three-level atoms with a cascade configuration, where the target state is a 
Rydberg state — a state with a large principal quantum number and thus large electric 
dipole moment — and the atoms are so close to each other that a Rydberg blockade [I] 
happens. 

Initially both atoms are in the ground state |0) and we have 1 00) = |0) £g> |0) for the 
initial state of the two-atom system. With the two driving lasers coupling in the same 
way and with the same strength to both atoms, the two-atom state will be invariant 
under the permutation of the atoms for all later times. We use a shorthand notation 
in which, for instance, 1 01) denotes the state with one atom in the ground state and 
the other atom in the target state, |01) = (|0) ® |1) + |1) ® \Q))/y/2. In total, then, 
six two-atom states participate in the evolution: 1 00) , |01), |0e), |le), |ee), and |11). 
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Referring to this order, the interaction-picture Hamiltonian is 



/ o 





Ht = h 






V o 



2 











71 













2A + 5 
A RB + 26/ 



(27) 



after the usual dipole and rotating-wave approximations have been applied. The energy 
of the double Rydberg state |11) is shifted by an energy ^Arb due to the strong dipole- 
dipole interaction between Rydberg levels. When |A RB | |f2 |,|fii|, the Rydberg 
blockade mechanism prevents both atoms from being excited to the target state at the 
same time, and this mechanism can be used to entangle atoms and implement quantum 
gates [SI E]- As a typical value of this energy shift, we take A RB to be a few times the 
size of A. Since the double Rydberg state |11) is not going to be populated, the goal is a 
population transfer from the double ground state 1 00) to the one-Rydberg-atom target 
state 1 01) , i.e., a transfer between the two states to which the two left columns and the 



two top rows in (27) refer. 



There are two pairs of resonant transitions: 1 00) -H- |01) and |0e) -H- |le). There 
are thus two different choices of relevant states, either only the two states 1 00) and |01), 
or the four states 1 00) , |01), |0e) and |le). Correspondingly, we have a choice between 
a 6 = 2 + 4 or a 6 = 4 + 2 split, as indicated in (27). It is also possible to apply a 



two-step Markov approximation, i.e. first approximate four states and in a second step 
approximate the final two states. But this has no additional benefit, it gives the same 
result as applying the Markov approximation directly to the 6 = 2 + 4 split. 

The plots in figure [3] show the differences between adiabatic elimination 
(zeroth-order Markov approximation with condition p4^,)) and first-order Markov 
approximation with condition ([24^,) as well as the impact of the choice of the number 
of relevant states. As one would expect, the approximation with only two relevant 
states (bottom plot) includes less information than the one with four relevant states 
(top plot). Indeed, in the case of two relevant states, the approximations give a smooth 
curve without the high-frequency components that one sees in the top plot because the 
states giving rise to these "wiggles" have been eliminated. The plots also confirm that 
the first-order Markov approximation gives reliable results for both choices of relevant 
states. 



5. Outlook 



Although we managed to give answers to all four of the questions asked in Section 2.2 
this is not the end of the story. It will be instructive to apply the first-order Markov 
approximation to other more complex situations where standard adiabatic elimination 
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Figure 3. Population probabilities for a two-atom Raman transition in a three- 
level cascade system with Rydberg blockade, as a function of At. The plots are for 
fio = 0.4A, Hi — 0.3A, 5 — and Arb = 5A, with four relevant and two auxiliary 
levels in the top plot, and two relevant and four auxiliary levels in the bottom plot. The 
curves that start near 1 show the probability of the double ground state |00); the curves 
that start at and increase to about 1 before decreasing again show the probability 
of the one-atom target state |01); and the curves that start near and never grow to 
large values show the probabilities of the other states. The black curves show the full 
solution, the dashed blue curves are for the adiabatic-elimination approximation, and 



the dotted red curves for the first-order Markov approximation with condition ( 24 1) 



is currently the method of choice. As a matter of general interest, the judicious selection 
of the interaction picture appears to be crucial and deserves further study. One should 
also extend the method to include spontaneous emission or, more generally, investigate 
how the method can be modified for open quantum systems. Lastly, we mention that 
a systematic approximation without the use of adiabatic elimination is possible as well 

0- 
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